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Abstract 

There are several (mathematical) reasons why Dupire's formula fails in the non-diffusion 
setting. And yet, in practice, ad-hoc preconditioning of the option data works reasonably well. 
In this note we attempt to explain why. In particular, we propose a regularization procedure of 
the option data so that Dupire's local vol diffusion process recreates the correct option prices, 
even in manifest presence of jumps. 

1 Failure of Dupire's formula in non-diffusion setting 

Local volatility models IHIZ], ''dS/S — aiodS, t)dW" , are a must-be for every equity option trading 
floor. A central ingredient is Dupire 's formula, which allovirs to obtain the diffusion coefficient ctioc 
directly from the market (or a more complicated reference model), 

<7lJK,T) = 2dTC/K^dKKC, (1.1) 

in terms of (call) option prices at various strikes and maturities0 Of course, there are ill-posedness 
issues how to compute derivatives when only given discrete (market) data; this inverse problem is 
usually solved by fitting market (option or implied vol) data via a smooth parametrization, from 
which (Jioc is then computed. 

On a more fundamental level, given an arbitrage-free option price surface 

{C{K,T) : > 0,T > 0}, 

there are two problems in a non-diffusion setting: 

• Lack of smoothness {OkkC is precisely the stock price density). For instance, the asymmetric 
Variance Gamma (AVG; cf. 4 ) process (Vj) has characteristic function 



,[exp(mVT)] 



l-ievu+{a'^v/2)u'^ 
const ■ u"^"^/" as u - 



T/u 
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^Throughout, we work under the appropriate forward measure to avoid drift terms. 
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When 2T/v > 1 (cquivalcntly: T > iy/2), the characteristic function is integrable on E and 
Xt admits a continuous density; the same is then true for St = Soex.p{XT + ^T). This is 
not the case, however, for T < v/2, and indeed the density (given exphcitly e.g. in [TTl p. 82]) 
has a singularity at the origin, so that C cannot possibly be twice continuously differentiable 
in K. 

• Blowup of the short end of the local volatility surface. Assume, for the sake of argument, that 

St = So exp{XT + ctWt + fJ-T), 

where X is a pure jump Levy process with Levy measure L — L{dx), and is a standard 
Brownian motion. For instance, in the Merton jump diffusion model, where Xt has normally 
distributed jumps, one can see 

T) w const • ^ SiS T 0, K ^ Sq. (1.2) 

In fact, this blowup can be fully quantified, thanks to a recently established saddle point 
formula |5| ; details are given in the appendix. There we show that the approximation (jl.2l) is 
also true in the NIG (normal inverse Gaussian) model, whereas in Kou's double exponential 
jump diffusion model the blowup is of order T~^/^. (The analysis is not restricted to expo- 
nential Levy models. Indeed, general jump diffusion resp. semimartingales are "tangent" , at 
a given point in space-time, to a Levy process, cf. [T31[T5] .) 

At least from a mathematical point of view, all this is strong evidence that p.ip is not meaningful 
in presence of jumps. In fact, if one views (jl.ip as a (forward) PDE for call option prices as function 
of K,T, then the analogous formula in a jump setting is a (forward) PIDE, which is thus the 
"natural generalization" of Dupire's formula in presence of jumps; cf. ij and the references therein, 
in particular towards so-called local Levy models [2]. 

For good or bad, practitioners use (jl.ip no matter what. Smoothness is usually not an issue 
since local vol is typically obtained from a smooth parametrization of market (implied vol) dataH 
On the other hand, the short-time blowup of local volatility, o'ioc(-,T) — ^ cx) as T — ^ 0, is an 
immediate obstacle, for it already makes it unclear to what extent there exists a (unique) strong 
solution to the stochastic differential equation dS/S = ctioc {S, t) dW, let alone how to sample from 
it. It is reported from practitioners that various ad- hoc truncation and mollification procedures 
are in place, say with parameter e, after which the local volatility is not explosive anymore; Monte 
Carlo simulations of ^^dS'^ /S'^ = af^^{S'^ ,t)dW'^ are then possible. There is every practical evidence 
that 

E[{S'T-K)+]^CiK,T) 

is a good approximation (for otherwise, risk management would not allow this in practice). And 
yet, to our knowledge, there has been no mathematical justification to date for this type of approx- 
imations. 

2 Regular izat ion of call prices: how to make Dupire work 

We now consider the situation of a given martingale {St) which creates a smooth call price surface, 
with (strict) absence of calendar and butterfly spreads. This situation is typical in the industry. 

^Doing so without introducing arbitrage is decisively non-trivial, see 1111 1121 . 
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2/rx^A_ 2dTCiK,T + e) 



(Of course, only option data is known, but upon a suitable parametrization thereof, one is in the 
just described situation.) In particular, ai^^{K,T) = 28x0 / K'^OkkC is well-defined, as long as 
T > 0, but may explode as T — >■ (thereby indicating the possibility of jumps). This situation is 
also typical in a generic non-degenerate jump diffusion settingH 

Theorem 1. Assume that (St) is a martingale (possibly with jumps) with associated smooth call 
price surface C , 

yK,T>0: C{K,T) =E[{St- K)^], 

such that OtC > and BkkC > 0, i.e. (strict) absence of calendar and butterfly spreads. Define 
e-shifted local volatility 

Then dS^/S"^ = a^iS"^ ,t)dW , started at randomized spot Sq with distribution 

P[S'q e dK]/dK = dKKC{K,e), 

admits a unique, non-explosive strong SDE solution such that 

Vif, T > : E[{SI^ - K)+] C{K, T) as e 0. 

Our assumptions encode that the model itself, i.e. the specification of the dynamics of S*, has 
regularization effects built in. The result is then, in essence, a variation of the arguments put 
forward in a recent revisit of Kellerer's theorem, see |14| and the references therein. 
Proof. By assumption, a1 is well-defined for all T > 0. (In general, i.e. without adding e and in 
presence of jumps, local vol is not well-defined in the sense that it may blow up as T — >■ 0.) Existence 
of a unique non-explosive local vol SDE solution for continuous and locally bounded diffusion 
coefficient is a classical result of the theory of (one-dimensional) SDEs. Set a^{S,t) = ^S^a'^{S,t); 
the generator of reads L ^ Lt ^ a^dss- Set also CiK, T) = C{K, T + e) and = C^ the 
density of Sr+e- By definition, 

^ ' ' p^{K,T) ' 
and hence, using dTKKC^{K,T) — 3tP^, 

dKK (a'p') = drp'. (2.1) 

In particular, p^ (which is (0, oo) n C°° (0, oo) in K) is a (classical) solution to the above Fokker- 
Planck equation. On the other hand, S'^ solves the martingale problem for Lt in the sense that for 
any test function ip, say, smooth with compact support, 

Jo 



^The mathematics here is well understood; a smooth density (and then call prices) can be the result of a 
(hypo)elliptic diffusion part, infinite activity jumps may also help. From a practical point of view, models tend 
to be locally elliptic with jumps super-imposed, so that we shall not pursue further technical conditions here. 
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is a zero-mean martingale. Taking expectations and writing q"^ — q'^{dS, t) for the law of noting 
g^(d5,0) =p^(5,0)dS', we see 

^{S)q'{dS,t)= J ^{S)q%dS,0) + J a'{S,s)^"{S)q'{dS,s), 

which is nothing but an (analytically weak, in space) formulation of the Fokker-Planck equa- 
tion p.ip . with solution given in terms of a family of measures, q'^ = {q^{dS,t) : t > 0}. By a 
suitable uniqueness result for such equations due to M. Pierre (see [H]), we see that p^dS = q'^ . 
And we have the following consequence for call prices based on e-regularized local vol: 

E[(5'|, - i^)+] = / {S - K)+ q''{dS,T) 

{S - K)+ p%dS,T) 



= mSr+e - K)+] 
= C{K,T + e). 

Since, by assumption, C is continuous, convergence as e — > to C (K, T) is trivial. ■ 

Remark 2. Instead of : T i-^ T + e one may take any strictly increasing T;, : [0,oo) [e,oo); 
in general, by the chain rule from calculus, an additional factor will appear in Dupire 's formula. 



3 Appendix: Local volatility blowup in some jump models 

Throughout this section, we normalize spot w.l.o.g. to 5*0 = 1. At the money, i.e. for K — Sq — 1, 
no blowup of the local volatility is to be expected. For instance, in Levy jump diffusion models, 
(Tioc tends to the volatility a of the jump diffusion part as T — >■ 0. To see this, recall the forward 
PIDE for the call price [H [21 H]: 

drC = -K^a^dKKC + / iy{dz) {C{Ke~', T) - C{K, T) - K{e' - 1)9kC) , 

^ J -oo 

where v is the Levy measure. As T — ^ 0, the integral tends to a non-negative constant. The claim 
thus follows from (|l.ip . since the density OkkC tends to infinity for K = Sq. 

By similar reasoning, we can quantify the off-the-money blowup of local vol, as soon as small- 
time asymptotics for the density in the denominator of Dupire's formula (jl.ip are available. For 
example, the density of the NIG (normal inverse Gaussian) model is ^ const ■ T (in fact, there is 
an explicit expression for the density 3 ). This implies a^^^{K,T) « const /T for K So- (Strictly 
speaking, the argument gives only a bound 0{1/T). To show that the numerator dxC tends to a 
nonzero constant seems difficult.) 

We now obtain further asymptotic results by using the local volatility approximation 

al,{K,T)^ 9^ \' / 3.1 
s(s — 1) 

^ ' s=s{K,T) 
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presented in [S1IH]- Here, 171(3, T) = log M(s,T) denotes the log of the moment generating function 
M(s, T) = i?[exp(sXT)] of the log-price Xt = St, and s = s{K, T) solves the saddle point equation 

d 

—m{s,T) = k:^ \ogK. 

While our focus in [5] was on asymptotics in K, may as well be used for time asymptotics, 

provided that the underlying saddle point approximation can be justified. The latter can be achieved 
by analysing the Fourier representation of local vol, which is at the base of 

'-loclA, ) - ^ ^^^^ . [6.2) 

For instance, in [8] we discussed the strike asymptotics for Kou's double exponential jump diffusion, 
and small-time asymptotics can be done in the same way, the result being that local variance is of 
order T~^/^. Note that tail integral estimates are required to give a rigorous proof; see [8j |9] for 
two examples of such arguments. For the Kou model, the moment generating function 

X+p , X-{l-p) 



M(s,r)=exp T^ + A 



Ai - s A_ + s 



has a singularity of the type "exponential of a pole" , and one can copy almost verbatim the saddle 
point analysis of On the other hand, the NIG model, which we discussed above, is not in 

the scope of (|3.ip . The reason is simply that there is no saddle point, as the moment generating 
function 



M{s,T) = exp 
has no blowup at its singularity a — j3. 

The Merton model 

The Merton jump diffusion is another case in point for the saddle point approximation (|3.ip . The 
resulting off-the-money blowup is essentially of order 1/T: 

\k\f, \k\^-^'^ 



aUK,T)^ const -^i^log^j , T ^ 0, (3.3) 

where k ~ logiiT 7^ and the constant depends on the sign of fc. To see this, recall that the Merton 
jump diffusion has moment generating function 

M(s, T) = exp{T{^ah^ + bs + \(e^"^"/2+f^s _ 

with diffusion volatility a, jump intensity A, mean jump size fi and jump size variance S^. The 
saddle point s solves dsm{s, 1) = k/T. As observed in p^;, for Levy models the saddle point (if it 
exists) is always a function oi k/T, so that various asymptotic regimes can be captured by the same 
formula. For fc > 0, the saddle point satisfies 



6^s^ ^ k V2^i L k 1^ , fc u2 ^ „/loglog 

= log — — \l log — - - log log — + -J - log — + O 
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as k/T — > oo. Inserting this into p.ip yields p.3p . with a similar reasoning for k <{). As regards 
the tail integral estimates needed to make this rigorous, no unpleasant surprises are to be expected, 
as double exponential singularities (here: at infinity) are well known to be amenable to the saddle 
point method. 



Merton's jump-to-ruin model 

In this model, the underlying [St) follows Black-Scholes dynamics, but may jump to zero (and stay 
there) at an independent exponential time with parameter A. In other words, A is the risk-neutral 
arrival rate of default. Out of the money, the additional default feature has little influence on local 
vol for small T, as it does not matter much that the underlying can jump even further out of the 
money. So we expect a'^^^{K,T) —)■ cr^ as T for if > 5*0 = 1. Indeed, this follows immediately 
from the explicit formula (see [lO l I16 j ) 

aL(if,T) = a^+2AaV^|M, (3.4) 

where N is the standard Gaussian cdf and 

, \og{So/K + \T) aVr 

Formula (|3.4I) is an easy consequence of the fact that the call price in this model is just the Black- 
Scholes price with interest rate A (A -f r in general, but recall that we assume r = throughout). 
The out-of-the- money convergence af^^{K,T) — > can also be confirmed by our formula p.l^ : 
The moment generating function equals 

M(., T) = + - ' ^ ° 

1 00 s < 0. 

There is a saddle point s = k/{a'^T) + \- A/ct^, and hence gives a^^^iK^T) = ct^ + o{T). 

For if < 5*0 = 1, on the other hand, (|3.4I) reveals a fast blowup of the order (Tj^Q^(if, T) w e^/"^ 
as T — ^ 0. Formula seems not to be useful here, as the moment generating function is not 

defined for Re(s) < 0, since St may assume the value zero. 
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